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ABSTRACT 

> 
00 . 

fSj . Using the multi-fluid approach, we investigate streaming and thermal instabilities of 

the electron-ion plasma with homogeneous cold cosmic rays drifting perpendicular to the 
background magnetic field. Perturbations across the magnetic field are considered. The 
back- reaction of cosmic rays resulting in the streaming instability is taken into account. 
The thermal instability is shown not to be subject to the action of cosmic rays in the 
model under consideration. The dispersion relation for the thermal instability has been 
derived which includes sound velocities of plasma and cosmic rays, Alfven and cosmic ray 
drift velocities. The relation between these parameters determines the kind of thermal 
instability from Parker's to Field's type instability. The results obtained can be useful for a 
more detailed the investigation of electron-ion astrophysical objects such as galaxy clusters 
including the dynamics of streaming cosmic rays. 
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1. INTRODUCTION 

The cosmic rays are an important ingredient in some of the astrophysical environments 
(e.g. Zweibel 2003). They are capable of affecting the dynamics of the astrophysical 
plasma media leading to plasma heating, increasing the level of ionization, driving outflows, 
modifying shocks, and so on (e.g. Field, Goldsmith & Habing 1969; Zweibel 2003; Guo 
& Oh 2008; Everett et al. 2008; Beresnyak, Jones & Lazarian 2009; Samui, Subramanian 
& Srianand 2010; EnBlin et al. 2011). Ionization by cosmic rays has a vital role in star 
formation near the Galactic center (e.g. Yusef-Zadeh, Wardle & Roy 2007) and in the dead 
zone of protoplanetary disks (Gammie 1996). 

The thermal instability (Field 1965) has been used to explain existence of the cold 
dense structures in the interstellar (e.g. Field 1965; Begelman & McKee 1990; Koyama 
& Inutsuka 2000; Hennebelle & Perault 2000; Sanchez-Salcedo, Vazquez-Semadeni & 
Gazol 2002; Vazquez-Semadeni et al. 2006; Fukue & Kamaya 2007; Inoue & Inutsuka 
2008; Shadmehri, Nejad-Asghar & Khesali 2010) and intracluster (ICM; e.g. Field 1965; 
Mathews & Bregman 1978; Balbus & Soker 1989; Loewenstein 1990; Bogdanovic et al. 
2009; Parrish, Quataert & Sharma 2009; Sharma, Parrish & Quataert 2010) media. For 
example, molecular filaments are seen in galaxy clusters with short (<1 Gyr) cooling times 
(e.g. Conselice, Gallagher & Wyse 2001; Salome et al. 2006; Cavagnolo et al. 2008; O'Dea 
et al. 2008). 

In galaxy clusters, cosmic rays are wide spread (e.g. Guo & Oh 2008; EnBlin et al. 2011). 
Therefore, they could exert influence on the thermal instability. In particular, including 
cosmic rays is required to explain the atomic and molecular lines observed in filaments in 
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clusters of galaxies (Ferland et al. 2009). Such an investigation has been performed by 
Sharma, Parrish & Quataert (2010) in the framework of the magnetohydrodynamic (MHD) 
equations where cosmic rays have been considered as a second fluid having the velocity of 
the thermal plasma. Numerical analysis has shown that the cosmic ray pressure can play an 
important role in the dynamics of cold filaments making them much more elongated along 
the magnetic field lines than the Field length. 

However, the relativistic cosmic rays can have a drift velocity of the order of the speed of 
light and temperature larger than the particle rest energy. The interaction of such particles 
with the thermal plasma can not be considered in the framework of the conventional 
MHD. The cosmic ray drift current results in arising of the return current provided by 
the background plasma (e.g. Achterberg 1983; Bell 2004, 2005; Riquelme & Spitkovsky 
2009, 2010). The possible role of this effect in the generation of thermal instability needs 
to be considered. There is also another important issue like the amplification of magnetic 
fields. The classical cyclotron resonant instability has been proposed long time ago to 
explain this process (Kulsrud & Pearce 1969). However, this mechanism has turned out to 
be unable to provide sufficient energy in the shock upstream plasma. In order to resolve 
this problem, just recently a new non-resonant instability has been proposed that may 
provide a much higher energy (Bell 2004; see also Zweibel 2003). This instability known 
as the Bell instability has been confirmed by non-linear numerical simulations (Riquelme 
& Spitkovsky 2009). Subsequent works extended this instability into various directions by 
considering partially ionized media (Reville et al. 2007) and thermal plasma effects (Zweibel 
& Everett 2010). However, the works sited above, except for the paper by Bell (2005), 
have been restricted to the cosmic ray drift velocity and perturbations parallel to the initial 
magnetic field. In his paper, Bell (2005) has derived the general dispersion relation for 
arbitrary orientation of the background magnetic field, cosmic ray current, and direction 
of perturbations. The dispersion relation obtained describes instability due to the return 
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plasma current. In the paper by Riquelme & Spitkovsky (2010), the case in which the 
cosmic ray current is perpendicular to the initial magnetic field has also been considered. 
In this case, cosmic rays can be magnetized in a way that their Larmor radius defined 
by the longitudinal thermal velocity (Zweibel 2003; Bell 2004) is smaller than the typical 
length scales of the system. Riquelme & Spitkovsky (2010) have studied this perpendicular 
current-driven instability analytically in the linear regime and numerically. Their growth 
rate was similar to that of the cosmic ray current-driven instability by Bell (2004). But 
these authors have not included the cosmic ray back-reaction analytically. 

The thermal instability in galaxy clusters in the multi-fluid approach has been 
considered by Nekrasov (2011, 2012). Effects related to cosmic rays were not included 
in these papers. Although the original Bell instability was proposed to explain magnetic 
field amplification in a shock, just recently Nekrasov and Shadmehri (2012) extended the 
instability to a multi-fluid case in which the thermal effects are also considered along 
with the presence of streaming cold cosmic rays. A geometry was considered in which 
homogeneous cosmic rays drift across the background magnetic field and perturbations 
arise along the latter. Such a geometry was analogous to that treated by Riquelme & 
Spitkovsky (2010). The cosmic ray back-reaction has been included and the growth rate 
has been obtained which is much larger than that for the Bell instability and perpendicular 
current-driven instability by Riquelme & Spitkovsky (2010). These findings motivated us 
to investigate the case in which perturbations arise transversely to the ambient magnetic 
field in the directions both along and across the cosmic ray current. As it is followed from 
the paper by Bell (2005) using the MHD equations, a streaming instability does not exist 
for such a geometry. However, this result is incorrect in the multi-fluid consideration that 
is shown in this paper and has been obtained earlier (see for example Nekrasov (2007)). 
We include the induced return current of the background plasma and back-reaction of 
cosmic rays. With such an approach, the dispersion relations are derived and the growth 



-5 - 



rates are found analytically. We also consider possible effects of cosmic rays on the thermal 
instability. For simplicity, we ignore the action of gravity as it has been done by Sharma, 
Parrish & Quataert (2010). The effects of the gravitational field have been investigated in 
detail by the multi-fluid approach in papers by Nekrasov & Shadmehri (2010, 2011). Thus, 
our present study extends previous analytical studies by considering not only the thermal 
effects but the currents driven by cosmic rays and their back-reaction. 

The paper is organized as follows. Section 2 contains the fundamental equations for 
plasma, cosmic rays, and electromagnetic fields used in this paper. Equilibrium state is 
discussed in Section 3. Wave equations are given in Section 4. In Sections 5 and 6, the 
dispersion relations including the plasma return current, cosmic ray back- reaction, and 
the terms describing the thermal instability are derived and their solutions are found 
for perturbations along and across the cosmic ray current, respectively. Discussion of 
important results obtained and possible astrophysical implications are provided in Section 
7. Conclusive remarks are summarized in Section 8. 



2. BASIC EQUATIONS FOR PLASMA AND COSMIC RAYS 

The fundamental equations for a plasma that we consider here are the following: 

^ + v, ■ Vv, = - + «LE+ JLv, x B, (1) 



dt m j n j m j m j c 



the equation of motion, 



the continuity equation, 



^ + V • n, Vj = 0, (2) 



dT 1 

_i + Wi . V T, + ( 7 - 1) TjV ■ v, = - ( 7 - 1) —d {m, Ti) + v\ e (n e , T e ) (T e - T t ) (3) 
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and 

dT 1 

+ v e • VT e + (7 - 1) T e V • Ve = - (7 - 1) ~Ce K, T e ) - ^ (m, T e ) (T e - T,) (4) 
at n e 

are the temperature equations for ions and electrons. In Equations (1) and (2), the index 
j = i,e denotes the ions and electrons, respectively. Notations in Equations (l)-(4) are the 
following: qj and rrij are the charge and mass of species j, Vj is the hydrodynamic velocity, 
rij is the number density, pj = UjTj is the thermal pressure, Tj is the temperature, vf e (n e , T e ) 
(v e ei (rii,T e )) is the frequency of the thermal energy exchange between ions (electrons) and 
electrons (ions) being uf e (n ei T e ) = 2u ie , where u ie is the collision frequency of ions with 
electrons (Braginskii 1965), riiuf e (n e ,T e ) = n e v e ei (rii,T e ), 7 is the ratio of the specific heats, 
E and B are the electric and magnetic fields, and c is the speed of light in vacuum. For 
simplicity, here we do not take into account collisions between the ions and electrons in the 
momentum equation. This effect for the thermal instability has been treated by Nekrasov 
(2011, 2012), where, in particular, conditions, under which such collisions can be neglected, 
have been found. However, the thermal exchange should be included because it must be 
compared with the dynamical time. The cooling and heating of plasma species in Equations 
(3) and (4) are described by the function Cj(rij,Tj) = n'jAj (Tj) — rijTj, where Aj and Tj 
are the cooling and heating functions, respectively. The form of this function has a certain 
deviation from the usually used cooling- heating function £ (Field 1965). Both functions are 
connected to each other via the equality Cj (n,j,Tj) = rrijnj£j. Our choice is analogous to 
those of Begelman & Zweibel (1994), Bogdanovic et al. (2009), Parrish, Quataert & Sharma 
(2009). The function A^- (Tj) can be found, for example, in Tozzi & Norman (2001). We do 
not take into account the transverse thermal fluxes in the temperature equations, which are 
small in the weekly collisional plasma (Braginskii 1965) being considered in this paper. 

Equations for relativistic cosmic rays which can be in general both protons and 
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electrons we use in the form (e.g. Lontano, Bulanov & Koga 2002) 

+ V cr • V (RcrPcr) = + qcr E+-V cr X B , (5) 



where 



oc n cr V c 

iZcr = 1 + ^-^~ 2 . (7) 



In these equations, p cr = 7 cr m cr v cr is the momentum of a cosmic ray particle having the rest 
mass m cr and velocity v cr , g cr is its charge, p cr = % r 1 n cr T cr is the kinetic pressure, n cr is the 
number density in the laboratory frame, T cr is the adiabatic index, 7 cr = (1 — v^ r /c 2 ) 
is the relativistic factor. The continuity equation is the same as Equation (2) for j = cr. 
Equation (7) can be used for both cold nonrelativistic, T cr m cr c 2 , and hot relativistic, 
T cr ^> m cr c 2 , cosmic rays. In the first (second) case, we have T cr = 5/3 (4/3) (Lontano, 
Bulanov & Koga 2002). The general form of the value R cr , which is valid for any relations 
between T cr and m cr c 2 , can be found e.g. in Toepfer (1971) and Dzhavakhishvili and 
Tsintsadze (1973). 

Equations (l)-(4), (5), and (6) are solved together with Maxwell's equations 

VxE = -i» (8) 



and 



„ „ 4tt. IdE 



where j = j pi + j cr = Ylj Qj n j y j + Jcr- Below, we first consider an equilibrium state in which 
there is a stationary cosmic ray current. 
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3. EQUILIBRIUM STATE 

We will consider a uniform plasma embedded in the uniform magnetic field B (subject 
here and below denotes background parameters) directed along the z-axis. We assume 
that the plasma in equilibrium is penetrated by a uniform beam of cosmic rays having 
the uniform streaming velocity u cr along the y-axis. It is reasonable to suppose that in 
such uniform model the magnetic field due to cosmic rays will be absent. This picture is 
analogous to the consideration of the gravitational potential in the equilibrium state in an 
infinite uniform medium having a constant mass density. Then we obtain from Equation (9) 

QjnjoVjo + j«-o+^^ = 0. (10) 

3 

From Equation (1), we easily find in the equilibrium state 

E x B c <9E 

v e0 = v = c — 2 — , v i0 = v H - — 7T-, v e0z = v i0z = 0, (11) 

-DO W« -DO at 

where we have assumed that d/dt <C u C j, u C j = qjB /mjC is the cyclotron frequency. 
Analogously, we obtain from Equation (5) under condition R cr ^f cr d/dt <C u cr (u cr = 
q cr B /m cr c) 

v cr0 = v + u cr . (12) 

In Equation (12), we have neglected the polarization drift of cosmic rays (the second term 
on the right-hand side for the ion velocity in Equation (11)). This can be done for the 
approximate condition n i0 ^> n cr o, if the cosmic rays are protons (see below), which is 
always satisfied. Using Equations (11) and (12), we will find the current j 

. _ qiniocdE 

JO — E~~i^r + Qcr n crO u cr, 

where we have taken into account the condition of quasi-neutrality qin iQ + q e n e0 + q cr n cr0 = 
(the number density n cr is the one in the laboratory frame). This condition is satisfied 
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in astrophysical plasmas due to cosmic ray charge neutralization from the background 
environment (Alfven 1939). Substituting Equation (13) into Equation (10), we obtain 

__ = -47rg cr n cr0 u cr 2 , (14) 
at c Ai + c 

where ca% = (Bq / 'inminio) 1 ^ 2 is the ion Alfven velocity Usually the inequality c 2 ^> c Ai is 
satisfied. In this case, Equation (14) coincides with the corresponding equation given in the 
paper by Riquelme & Spitkovsky (2010). Substituting Equation (14) into Equation (11) for 
ions, we find the return plasma current 

iret = Qi^io (Vjo — v ) = — q cr n cr0 u cr , (15) 

which magnitude is equal the cosmic ray current and has the opposite direction. From 
Equation (15), it is followed that the induced plasma current drift velocity u p i = v i0 — v 
is equal to u p ; = — (? C r^cro/%^io) u cr- Using Equation (14), we see that the polarization 
cosmic ray drift velocity (-R cr 7 cr c/cu cr -B ) dE /dt can be neglected in comparison with u cr 
under condition 3> R cr ^ cr m cr n cr Q. If the cosmic rays (protons) are not too relativistic, 

this condition is satisfied. Below, the plasma drift velocity u p i will be also taken into 
account as u cr . 

We will consider the case in which background temperatures of the electrons and ions 
are equal, i.e. T e0 = T i0 = T . However, to follow the symmetric contribution of the ions and 
electrons in a convenient way, we make the calculations by assuming different temperatures. 
In this case, the thermal equations in equilibrium are given by 



A (nio, Tjo) — C e (n e0 , T e0 ) — 0. 



(16) 
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4. WAVE EQUATIONS 

For perturbations across the background magnetic field, Equations (8) and (9) give us 
the following two equations: 



and 



where ji = j p a + j cr i and the subscript 1 here and below denotes the perturbed values. 
The third equation describes the ordinary electromagnetic wave with Ei|| B Q . The general 
expressions for the components j p ii Xty and j cr ix, y are given in the Appendices A and 
B (Equations (A54)-(A56) and (B19)-(B21)). These expressions are available for both 
magnetized and non-magnetic systems, electron-positron, pair-ion, and dusty plasmas and 
so on. In their general form, they are very complicated. Therefore to proceed analytically, 
one must apply simplifying assumptions. We are interested in magnetized systems consisting 
of electrons, ions, and cosmic rays where cyclotron frequencies of species are much larger 
than inverse dynamical times. In our case, this implies 

^(l+^l) 2 ' (19) 

2 4 fd d\ 2 

(see Equations (A5), (A8), and (B7)). As we have noted above, the cosmic rays can be both 
protons and electrons. For ultrarelativistic cosmic rays, 7 cr0 1, the second inequality 
(19) can be violated. Such a case is not considered here. Another condition that simplifies 
the treatment considerably is to assume the wavelength of perturbations to be much larger 
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than the Larmor radius of particles pj 

1 » p 2 V 2 , (20) 
2 / 2 d 2 d 2 



dx 2 dy 2 / 

(see Equations (A41) and (Bll)). The additional inequalities for cosmic rays will be 
given below. The third simplification is to consider perturbations along and across the 
cosmic ray drift velocity separately. The first case is simpler. Therefore, we begin with its 
consideration. 



5. THE CASE ^ ^ 0, £ = 

Using Equation (A56) and performing calculations of the corresponding quantities, we 
find the components of the plasma dielectric permeability tensor (v i0y has been changed by 

Upl) 



£plxx c4 \dt +Upl dy ' 17 

,2 



Eplxy 



Si J_ 

u 2 ci mi 



UJ 2 iUJ c 



Tin + T e0 — 



d_ 

gi + G 3 d 
D dt 



G 2 + G 4 ( d 



D 



u% mi 



+ 



Tic 



G 2 + G 4 ( 8 



D 



d 



£plyx 



a 2 



u 3 ci mi 



-i0 



D dt 



d_ 

dt 

D 



d 



dt +Upl dy 



d 



d 



dt +Upl dy 



dy 2 



£plyy 



2 

pi 

J: 

CI 



(21) 



df [dt 
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In obtaining expressions (21), we have taken into account that rrii ^> m e and n i0 ~ n e o- 
Analogously from Equation (B21), we obtain the cosmic ray dielectric permeability tensor 



2 / o o \ 2 / o \ -2 2 

per 3 / l / l per 2 2 



£cra:C ~ ^ 7cr ° V«9t + Ucr fy J J ~ ul r lcr ° Cscr I - 2 + ) a.. ( J > ( 22 ) 



■u cr (9 








c 2 ^ 




1 dy 1 





_ U pcr f d_\ U pcr 3 2 / Mcr g g \ g / g 

^ - ecrwx - ^ w «r ^ J + w 3 cr 7cr ° Cscr V c 2 dt dy J dy \dt 



cr \ / ccr 

2 



_ ^pcr 
£-cryy 2 IcrO- 

^ccr 

Here, we have used the additional condition for cosmic rays 

, 2 u cr ( d d \ d 

(see Equation (Bll)). The contribution of the term proportional to u cr /c 2 in Equation (22) 
is small. This term shows the contribution of 7 cr i to the cosmic ray pressure perturbation 
(see Equations (B8) and (B9)). 



5.1. Wave equation 

From Equations (17) and (18), using Equations (A54), (A55), (B19), and (B20) and by 
omitting the contribution of the displacement current, we obtain the equation 

2 fd\~ 2 d 2 E lx _ ( 

£ yy C I J Qy2 \ e xx£yy £ xy £ yx) ^\xi \ Z "->) 

where = e pli j + e cri j. The values are defined by Equations (21) and (22). When 
calculating the right-hand side of Equation (23), we assume some additional conditions 
except those given by Equations (19) and (20). We will neglect the contribution to e xy e yx 
of the thermal cosmic ray term in e crxy and e cryx . Besides, we will use the condition of 
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quasineutrality in e xy and e yx and neglect the terms arising due to expansion of fi i cr . An 
analysis shows that the corresponding conditions can be written in the form 



® d \ 2 2 ( u cr d t d \ d 



max < 7 -° U + Ucr d~y ) ; ° scr Wdi + d-y)d- y < (24) 

3 c\ cr f u cr d d \ 2 d 2 c 2 spl c 2 scr fu cr d d \ <9 3 



" lcr0 u 2 ccr \c 2 dt + dy) dy 2 ' lcr{) uj cl uj ccr \c 2 dt + dy ) dy*> 
c 2 ( u„ d d \ d ( d d x " 



lcr0 u ci uj ccr ^ c 2 dt + dy ) dy \ y dt + Upl dy 

1 II 

where c sp i = (2^Tio/mi) . According to conditions (20) and (24), the contribution of the 
term e xy e yx to the Equation (23) is small. Thus, we obtain 

= e xx ( E lx . (25) 



dy 2 xx \dt 



5.2. Dispersion relation 



Using Equations (21) and (22) and accomplishing the Fourier transform in Equation 
(25), we find for perturbations of the form exp (ik y y — iut) the following dispersion relation: 



,,2 , ,2 

Lt/ ™ - -2 , M pcr 3 / , n2 

IcrO (W - k y U cr ) 



= —(u- k y Upl) + 
,2 



(26) 



U pi h 2 1 

u 2 ci y mi 



Tio + T e0 H iu H % (u - kyUpi) 



D 



D 



UJ. 



per 2 7,2 2 r,'Z Z 

1crO li y C scr ~ ^y C • 



2„2 



Below, we consider solutions of Equation (26) for the streaming instability and an influence 
of the streaming and thermal pressure effects on the thermal instability. 
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5.2.1. Streaming instability 

Let us set all frequencies Vl equal to zero in Equation (26). To be more specific, it 
means that uj — k y u p i ^> Qt,™, and u ^> fi rne , Q t , where Q ie ~ fi ei = fi e (the frequencies 
f2 are defined by Equation (A12)). Then, this equation takes the form 

= (u - k y u p i) 2 + ^7cro ( w - K u crf (27) 

2 2 \ 

Upi 9 , V 2 2 , 2 \ >2 

The solution of Equation (27) is the following: 

" = ^ ± ^ [- - I*) 2 d+(l + d) (0%, + 7c ">L + 4)] V2 , (28) 

where 

2 2 

« — 2 2 y cr0 — y cr o- 

We see that the streaming instability has a threshold defined by the sound and ion 
Alfven velocities 

Ulrth ={l + d~ l ) (C 2 spl + 7c~rO rfc L + 4i) ■ (30) 

When this threshold is exceeded, u 2 r ^> u 2 rth , the growth rate 5 gr is given by 



d 1 ' 2 

$gr = ^ ^kyU C r- (31) 

These perturbations move with the phase velocity = (u p / + d-u cr ) / (1 + d). 



5.2.2. Thermal instability 



We now take into account the terms describing the thermal instability in Equation 
(26). We consider the fast thermal energy exchange regime in which fl £ 3> d/dt, &Ti,e, 
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Using Equations (A29) and (A30), we will have 

7 (2u - kyU p i) + in T , a _ ^_ 2 ( Aj2 __i^ 2 ^2 , ^ 2 
7 (2u - kyU p i) + i^VL T 

where 



C spl (d U lr - 7crO dc L - CM + ^ , (32) 



When obtaining Equation (32), we have assumed that uj <C k y u cr . If the right-hand 
side of Equation (32) is much less than unity, we obtain Field's isobaric solution 
2u = kyU p i — i£tr,n/l (Field 1965). These perturbations travel with the phase velocity 
Upi/2. In the opposite case, Equation (32) has Parker's isochoric solution 2u = k y u p i — iflr 
(Parker 1953). Thus, the presence of streaming cosmic rays can change the kind of thermal 
instability. When the right-hand side of Equation (32) is of the order of unity, the limiting 
solutions intermix. 



- 16 - 



6. THE CASE f- ^ 0, f- = 



Calculating the components of the plasma dielectric permeability tensor given by 
Equation (A56), we obtain 



-plxx 



-pixy 



UJ 
UJ. 



2 

pi 

2~' 

ci 



2 



+ 



UJ. 



u> 



n 2 

,2. 

pi 



+ 



J 2 t 

pe 



1 

rrii 



£plyx 



2 

Upi^ci 



a 2 



Tm — 

d_ 

dt 



G 3 + G 4 d 



D dt j dx 



d 





fd\ 


— 

dx 


Kdt) 



e piyy 



UJ. 
UJ. 



2 
pi 

3~ 

ci L 
2 

pi 

T 

ci 
2 



— IT, 



mi 



G 4 d 



■iO 



9_ 

ft 

Gi + G2 + G3 







— 

dx 


(1) 



. Si 2 <9 2 
H — -u r 



d 



Tio + T e0 — 

-2 



,2 

pi 

ci 



D 

1 G 2 -G 3 <9 3 



d 



dt J dx 2 V dt 



uj 2 ci pl dx 2 \dtj u^rm D d.r>\dt 
From Equation (B21) for cosmic rays, we will have 



-crxy 



-cryx 



-cryy 



t ; 2 

iu pcr 3 

2 icrOi 
ccr 








dt) 


w cr V 




UJp Cr UJ ccr 


(dt 


tt 2 

cr 





-1 



UJ, 



per 3 



9 



7cr0 I 7cr0 C scr a ^ccrU C r ] a I a , 



9 /9 



per 3 



UJ 



3 TcrO I ^scr 



d_ 

dx 



"t~ UJ ccr U c 



dx V eft 



9 /<9 



UJ 



per 



UJi 



IcrO 



1 + 7, 



d 2 (d 



crO^cr 



dx 2 V dt 



2 r*2 

^pcr 2 2^ 
- '-'scr^'crO 



d 



UJ 



ccr 



dx 2 \ dt 



-2 



(33) 



(34) 



In this case, the additional condition for cosmic rays except for Equation (20) is the 
following: 



1 >7, 



t^cr^s 
crO ^2~ 



d_ 
~ pcr dx~ 



(see Equation (Bll)). 
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6.1. Wave equation 

In the case under consideration, the wave equation has the form 

&xxC ^Qf^ Ox 2 i^xx^yy &xy£yx) -^ly (^5) 

Using Equations (33) and (34) and calculating the right-hand side of Equation (35), we find 

( - \- f i 4 1 , i 4cr\_ i 4c L 2 2 9 2 f ay 2 

\^xx e yy £xy £ yx) £ ra I y 2 ij 1 7cr0 J £xx ^2 ^crO C scr \ Qf I ^ ' 

<±( T . MT G 1 + G 2 + G 3 + G 4 d\ d 2 ( d^- 2 
^ulmX' ° D dt)dx 2 \dt, 

2 ,,,2 Q2 / Q\ -2 



. "pi "per 3 / n2 9 2 ( 8 \ 



6.2. Dispersion relation 

After Fourier transformation of Equation (35) and substitution of Equation (36), we 
derive the dispersion relation 

2 2 \ 2 2 

UJ ■ UJ I UJ ■ UJ 

pi ^ccr / pi ^ccr 

..a i f T , ^1 + ^2 + ^3 + ^4 . 
+ fc^ — I Tj + T e0 H 



III 

,2 



1 U> pcr 3,2/ _ \2 
lcrO K x \ u ct Upl) ■ 



£xx ^ccr 



Below, as above, we consider the streaming instability and influence of cosmic rays on the 
thermal instability. 
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6.2.1. Streaming instability 



As above, we again neglect in the values G^i = 1,2,3,4, and D all the frequencies Q 
Then, Equation (37) takes the form 

' 2 d . 



(1 + Icrld) 



C spl + lcrQ^ C scr + C Av 



(38) 



[l + d) ^^^spl 

where we have omitted u p i in comparison with u cr . This equation describes an aperiodic 
instability, if the drift velocity of cosmic rays exceeds the threshold value given by Equation 
(30). An expression for the growth rate 5 gr when u cr exceeds u crt h is the following: 

1/2 

kxUcr- (39) 



Jgr 



d 



(1 + d) (1 + 7-qgO 



6.2.2. Thermal instability 



Now, we take into account the contribution into Equation (37) of terms describing the 
thermal instability in the fast thermal energy exchange regime Q t ^> d/dt,Q,Ti,e,^ni,e- The 
dispersion relation becomes 



27a; + iQ>T,i 



d 



2 7 w + i^VL T spl 

This equation is analogous to Equation (32). Depending on whether the right-hand side 
of Equation (40) is much larger than the unity or not, we will have Parker's or Field's 
instability (see above). In these limiting cases, the value u 2 on the right-hand side of 
Equation (40) must be substituted by — Vl\ or — f^„, respectively. 



(40) 
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7. DISCUSSION AND IMPLICATIONS 

The growth rates (31) and (39) of streaming instabilities have a similar form and 
increase with decreasing of the perturbation wavelength. The thresholds for the cases 
k x = 0, k y 7^ and k x ^ 0, k y — are equal to each other (see Equations (28) at u cr ^> u v \ 
and (38)). Thus, streaming cosmic rays generate perturbations in all directions across 
the ambient magnetic field. A spectrum of the perturbations in the k-space is limited 
from above by conditions given by Equations (19) and (20) and additional conditions (see 
inequalities after Equations (22) and (34)). These conditions for the case k x — 0, ky ^ can 
be written in the form 



\ y \ j d u cr _ 7 cr0 u cr _ 7 cr0 u cr c 

^> max 



where the value d is defined by Equation (29) and A is the wavelength. We have assumed 
that the threshold of instability is exceeded. The conditions (20) are satisfied. The 
analogous conditions for the case k x ^ 0, k y = are the following: 



^ ma¥ j ^ ^cr. ^crpd ulr . 2 3 2 4 u lr <4cr 

1^ I s , — 2j\ , ,2 ' /i , j\ /i , -2 j\ , ,9 ' "» ' IcrOFcri IcrO, ,o „A 

w ccr ccr 

Let us consider the polarization of perturbations. In the case k x = 0, k y ^ 0, the current 
jiy = (see Equation (18) without the displacement current). Then, the component of the 
electric field E iy is equal to E ly = — (s yx /e yy ) E ix . Estimations show that (s yx /e yy ) «C 1 
for the streaming instability. Thus, the polarization is a linear one being E ix 3> E ly . 
The electric field polarization for the thermal instability depends on the wavelength of 
perturbations, ion Alfven velocity, and parameters of cosmic rays and can be various. In the 
case k x ^ 0, k y = 0, the current j lx = (see Equation (17)). Then E lx = — (e xy /e xx ) E ly . 
The ratio e xy /e xx for the streaming instability is given by 

1/2 

+ terms <C 1, 



d{l + l- 2 d) 



(1 + d) 
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where f — 1 — qcrmi/lcroQi 171 ^- We have assumed that 

j 2 kx C scr ^ j 

(1 + d) crC W r <5 3r 

where <5 flr is determined by Equation (39). This condition can be easily satisfied taking into 
account Equation (20). Thus, we see that the ratio e xy /e xx can be both smaller and larger 
than the unity. This result is also just for the case of thermal instability. 

From Equations (32) and (40), it is followed that the relations between hydrodynamical 
parameters of thermal plasma and cosmic rays and the perturbation wavelength determine 
the kind of thermal instability from Parker's(1953) to Field's (1965) type instability. 
Dissipative processes such as the thermal conductivity of plasma and cosmic rays can affect 
the growth rate of thermal instability. 

We now compare the growth rate found for the streaming instability along the 
background magnetic field (Nekrasov & Shadmehri 2012) with the growth rates obtained 
in this paper. The growth rates given by Equations (31) and (39) are of the same order 
of magnitude, if 7 cr o ~ 1 or 7 cr0 1 and d < 1 (for the same wavenumbers). In the case 
Icro 3> 1 and d 3> 1, the growth rate given by Equation (39) is larger. Therefore, we use 
Equation (39) for a comparison. The maximal growth rate found by Nekrasov & Shadmehri 
(2012) is equal to 

r _ ■ I n \ ( lcrO C A 



m cr n cr0 c 2 J \l cr l c 2 cr . + c 2 A 
where ca — ca% (l + %ro^) ■ The ratio of this growth rate to the growth rate (39) for the 
same cosmic ray drift velocities is the following: 

&m /-, . ,_i\l/2 CAi ^pcr 

— = {i + d 1 ) mir~ ■ 

V (c 2 cr + -f cr0 c 2 A ) k xC 

We see that for sufficiently short wavelengths the ratio S m /S gr can be less than unity. Thus, 
the transverse streaming instabilities induced by cosmic rays can considerably contribute to 
turbulence of astrophysical objects and amplification of magnetic fields. 
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We have explored the situation in which cosmic rays drift across the background 
magnetic field. This model has been considered by Riquelme & Spitkovsky (2010) for the 
problem of the magnetic field amplification in the upstream region of the supernova remnant 
shocks. The perturbations along the background magnetic field have been investigated and 
the cosmic ray back- reaction has not been taken into account in the analytical treatment. 
The latter effect for the longitudinal perturbations has been included in the paper by 
Nekrasov & Shadmehri (2012) where the growth rate considerably larger than that in 
(Riquelme & Spitkovsky 2010) has been found. In this paper, we have investigated the 
transverse perturbations. In the paper by Bell (2005), the unstable perturbations for the 
last case are absent in the MHD model. However, the multi-fluid approach gives a different 
result. In another model, cosmic rays drift along the magnetic field. This case has been 
investigated by Bell (2004) (see also Riquelme & Spitkovsky 2009). In both cases (Bell 
2004; Riquelme & Spitkovsky 2010), the growth rates have turned out to be the same. 

The streaming cosmic ray driven instabilities can exist in a variety of environments. 
Although such a type of instability was suggested originally for the magnetic field 
amplification in the shocks of supernovae, we think, wherever there is a strong cosmic 
ray streaming, this instability may play a significant role. For example, the models 
described above can be applied to the ICM where cosmic rays are an important ingredient 
(Loewenstein et al. 1991; Guo & Oh 2008; Sharma et al. 2009; Sharma, Parrish & Quataert 
2010). Observations show that many cavities or bubbles in the ICM contain cosmic rays and 
magnetic field (e.g. Guo & Oh 2008). A substantial amount of cosmic rays may escape from 
these buoyantly rising bubbles (e.g. EnBlin 2003) which could be shredded or disrupted by 
RT and KH instabilities as they rise through the ICM (e.g. Fabian et al. 2006). Cosmic 
rays may also be produced by other processes near the central AGN of the galaxy cluster. 
Structure formation shocks, merger shocks and supernovae may also inject cosmic rays into 
the ICM (e.g. Voelk, Aharonian & Breitschwerdt 1996; Berezinsky, Blasi & Ptuskin 1997). 
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The observation of diffuse radio synchrotron emission in many galaxy clusters give direct 
evidence for the presence of an extensive population of non-thermal particles (e.g. Brunetti 
et al. 2001; Pfrommer & EnBlin 2004). Recent Chandra and XMM observations also show 
evidence for a significant non-thermal particle population within the ICM (Sanders, Fabian 
& Dunn 2005; Werner et al. 2007). 

In some of supernova remnants such as IC 443, SN 1006, Kepler, Tycho and etc., 
the driven shocks are propagating in a partially ionized ambient medium. This was a 
good motivation to extend cosmic ray streaming instability from the MHD approach to 
a two-fluid case, by considering ions and neutrals as two separate fluids where they can 
exchange momentum via collisions (e.g. Reville et al. 2007; see also Bykov & Toptygin 
2005). It has been shown that the instability is getting slower rate because of collisions 
of ions with neutrals, in particular when the cosmic ray flux is not very strong. However, 
the back-reaction of cosmic rays has not been considered. Having in mind the finding 
that the growth rate is significantly enhanced in the presence of cosmic ray back-reaction 
in a three-fluid plasma system consisting of the ions, electrons, and cosmic rays, one 
may naturally expect such an effect in a four-fluid plasma system consisting of the ions, 
electrons, cosmic rays, and neutrals. It deserves a further study, but we may expect that 
the stabilizing effect of the ion-neutral collisions can be compensated by the back-reaction 
of cosmic rays. 



8. CONCLUSION 

Using the multi-fluid approach, we have investigated streaming and thermal instabilities 
of the electron-ion plasma with homogeneous cold cosmic rays drifting across the background 
magnetic field. We have taken into account the return current of the background plasma 
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and the back-reaction of cosmic rays for perturbations transverse to the magnetic field 
and along and across to the cosmic ray drift velocity. For sufficiently short wavelength 
perturbations, the growth rates exceed the one of streaming instability along the magnetic 
field. 

The thermal instability has been shown not to be subjected to the action of cosmic 
rays in the model under consideration. The dispersion relations for the thermal instability 
in the multi-fluid approach has been derived which include sound velocities of plasma and 
cosmic rays, Alfven and cosmic ray drift velocities. The relations between these parameters 
determine the kind of thermal instability from Parker's to Field's type instability. 

The results of this paper can be useful for the investigation of the electron-ion 
astrophysical objects such as galaxy clusters including the dynamics of streaming cosmic 
rays. 
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A. Appendix 

A.l. Perturbed velocities of ions and electrons 

We put in Equation (1) v,- = v j0 + v^i, Pj = Pjo + Pji, E = E + Ei, B = B + Bi, 
where the subscript denotes equilibrium uniform parameters and the subscript 1 relates 
to perturbations. Then the linearized version of this equation takes the form 

-7— + Vj • VVji — h r jiH Vji X ti , (A1J 

at m,j rrij rij rrijC 

where we have used that pji = rijoTji + n^Tjo {rij = n j0 + nji, Tj = T j0 + Tji) and 
introduced notation 

F = iL El + Jz- V x Bx. (A2) 
rrij rrijC 

From Equation (Al), we find expressions for the ion velocities Vn Xjy in the form 

1 T rr T t ^ fd d\- 1 



^iVn x = — L ix Tn L ix — + v i0y — V • Vii (A3) 

rrii rrii \dt ay / 

fd d \ 

+U c iFny + I — + v iOyQ^ I *zlx 



and 



^ = ^ l ^-^ l -(I + ^4) v ' v " (A4) 

In Equations (A3) and (A4), we have used the linearized continuity equation (2). The 
following notations are here introduced: 

n H^°4) 2+u4 (A5) 

d fd d\ d 

Li X — —UdT; 777 + Vio y — — , 



<9y \ <9t dy J dx 



d fd d\ d 

Jiy - Uci fa-\di + Vi0y dy) d~y 
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Analogous equations for the electrons are the following: 

tt 2 e v elx = — L ex T el - — L ex ( J V-v el +w ce F ely + ^^, (A6) 
m e m e \ot J at 

1 T" o f d \ OF i 

&e V ely = —L ey T e l — —L ey i— J V ■ V ei — UJ ce F elx -\ 7^—, (A7) 



where 



^ = ^ + b&, (A8) 

__ d__ JP_ 



9 <9 2 



Ley u ce ^ dydt 



We do not consider the longitudinal velocity Vji z because as can be shown in the case 
d/dz = this velocity only depends on the electric field E lz , dvj lz /dt = (qj/mj) E lz , and 
the transverse and longitudinal wave equations are split. 



A. 2. Perturbed temperatures of ions and electrons 

We find now equations for the temperature perturbations T i el . We here assume 
that equilibrium temperatures T i0 and T e0 are equal one another, T i0 = T e0 = T . The 
case T i0 ^ T e0 for thermal instability has been considered by Nekrasov (2011, 2012). 
For equal temperatures, the terms connected with the perturbation of thermal energy 
exchange frequency in Equations (3) and (4) will be absent. However for convenience of 
calculations, we formally retain different notations for the ion and electron temperatures. 
From Equations (3) and (4) in the linear form, we obtain equations for the temperature 
perturbations 

D u Tn - D 2l T el = C u V • Vi i , ( A9) 

D le T el - D 2e T a = C le V • v e i, (A10) 
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where notations are introduced 
D u = 



9 d \ o 

fa + V iOVQ- J + ^Ti + Ik 

fd d\ 



d 



d 



(All) 



Cu — Ti 



/9 \ 9 

D 2e = VL e i — i 



C\ e — T e o 



(7-D-|+a 



When obtaining Equations (A9) and (A10), we have used Equations (2) and (16). The 
frequencies in Equation (All) are the following: 

»Tj = (7 - 1) — — , Mnj = (7 - 1) — ^r^i , ( A12 ) 



n j0 dT j0 



T j0 dn j0 



^ie — v t e (^eOj 2~k)) , ^ei — ^li (^iO, ^eo) • 



From Equations (A9) and (A10), we find equations for Tn and T e i 



DTii = G 4 V • v a + G 3 V • v el 



(A13) 



and 



Here, we have 



DTel = GlV-Vei + G 2 V-Vii. 



(A14) 



L> = D u Di e - D 2i D 2e , 
G\ = DuC\ e ,G 2 = D 2e Cu, 
G3 = D 2 iC\ e , G4 = D\ e C\i. 



(A15) 
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A. 3. Expressions for V • v^, 



el 



We now substitute temperature perturbations T i>e i denned by Equations (A13) and 
(A14) into Equations (A3) and (A4). Then applying operators d/dx and d/dy to Equations 
(A3) and (A4), respectively, and adding them, we find equation for V • Vji 



LuV ■ v a = -L 2i V ■ v e i + 



where 



L>2i 



a 2 + - 



rrii 



~D [di + Vl0y dy~ ] ~ Xl ° 



V 2 , 



dFny dFn x 











dx dy J ' \dt Vl0y dy 



Analogously, using Equations (A6) and (A7), we obtain 

L le V • V e i = ~L 2e V ■ Vii + $ e l, 



where 



- T e0 V 2 , 



J 2e 



m P D dt 



. dF elv dF elx \ d ,_, ^ 



9a; cfy / dt 
From Equations (A16) and (A18), we find 



and 



The operator L is given by 



LV • v n = L u $n - L 2i <$> el 



IN ■ V e i = Lii$ e i -L 2e $ a . 



(A16) 



(A17) 



(A18) 



(A19) 



(A20) 



(A21) 



L — L\iL\ e — I 2 il 2e . 



(A22) 
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A. 4. Equations for ion and electron velocities via F 



i,el 



Using Equations (A3), (A4), (A13), (A20), and (A21), we obtain the following 
equations for components of the perturbed ion velocity: 



n2 



and 



ilx 



ITLiDL 



■>iy 



ITLiDL 

The operators Ai j2i are given by 
A ±i = 



I2i 



f d d\ 

(A^a - A 2i $ el ) + u cl F lly + f — + v l0y —\ F ilx (A23) 

f d d\ 

(A u $n - A 2i $ el ) - u ci F ilx + f — + Vio y —\ Fa v . (A24) 



-iO 



- iO 



d 


d 


dt 


+ v i0y - 


d 


d 


di 


+ v *»djj 



L 2 i — G^Lu. 



(A25) 



Equations for components of the perturbed electron velocity are found by using 
Equations (A6), (A7), (A14), (A20), and (A21) 



$£v elx 



Lrv Mi.'I',: - A 2e <f> n ) + u ce F ely 



m e DL 



dt 



L OF 
^Veiy = ( A ie$ei - A 2e <Z> n ) - 0J ce F elx + ' Ul 



m e DL 



dt 



Here, 



Ai e = 



A 2e = 



eO 



dt 



Lu — G 2 L 2i , 



L 2e — G 2 L\ e . 



(A26) 
(A27) 



(A28) 
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A. 5. Expressions for D and Gi, 2,3,4 

We now give expressions for D and Gi,2,3,4 denned by Equation (A15). Using Equation 
(All), we find 



¥t +v ^) U 1 D = 



dt +V * 0y &y ]+QTi 



9 o 

at 



(A29) 



9 



( — + fi Te ) fi ie 4 



9 d\ ' 



and 



eO 



9 9 

+ v iOyQ- ) + ^Ti + ^ie 



9* 



-(7-i)|+a 



9 9 



dt 



G2 — Tio£l ei 



C3 — TeO^je 



5 



-(7- 1) ( ) +^ni 



d_ 



9 9 



at 



— Tj ( — + fire + ^ei 



9f 



9 d 



9*' 



(A30) 



A. 6. Simplified expressions for A lj2 i and A 12e 



We can further simplify expressions for A\^i and A\^ e given by Equations (A25) and 
(A28). Using Equation (A17), we obtain 

(A31) 



A 2l = -G 3 Q 2 t . 



The expression for An can be given in the form 











G4 - DT i0 ( — + v i0y — 



-1 



n 2 e - —v 2 k, 



dt 1 ~™ y dy 

where we have used Equation (A19). The following notation is introduced in Equation 
(A32): 



(A32) 



K = I (C 2 C, - Gl G t ) (| + " | 



(9 



9 



-DT i0 T e0 [ — + v i0y ^ 



(A33) 
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Analogously, we will have 



A 2e = -G 2 n 



and 



mi 



d 



d \ f d 



W-—V 2 - + v Mv —\ - 



dt ' " i0y d y J \dt 



K. 



(A34) 



(A35) 



Calculations show that the value D 1 (G2G3 — G1G4) takes the simple form 



— (G2C3 — G1G4) — — T i0 T e0 



( d d\ 
(7-1) {ot+vvyQjj+nr- 



(A36) 



Using Equations (A29), (A30), and (A36), we can also write the value K defined by 
Equation (A33) in the simple form 



K = -T i0 T e0 (WiW e + WiQ ei + W e Q ie ) ^. 



(A37) 



Here, notations are introduced 



(A38) 



( d d\ 
Wi = 1 \ Qt + Vi0y dy J + ~ 
d 

W e = l^+ ^Te ~ fine- 

We remind the reader that the temperatures of the ions and electrons are considered to be 
equal one another. We retain different notations for the control of the symmetry of the ion 
and electron contribution. Analogously, we find the following values: 











G 4 - DT i0 ( — + v i0y — 



dt 



G 1 - DT ( 



eO 



dy 
d_ 
dt 



-1 



(A39) 



d 







where 



-T e0 (W e Vi + W e a ie + ViQei) [ — + v i0y — 



d d 



(A40) 



d_ 
di 



v e = — + n Te . 



A. 7. Operator L 



Let us find the operator L given by Equation (A22). Using Equations (A17) and (A19), 
we obtain 



1 



d 



G 4 ( d d 



— Tj 



iO 



n2 1 ^ 2 (G x d 

V H fr — — 

m P l \Ddt 



- T e0 1 V 



(A41) 



TYljUlpD V dt 



ViOy 



d_ 

dy 



V 4 K. 



The expressions containing in this equation are given by Equations (A37)-(A40). 



A. 8. Simplified equations for ion and electron velocities via Ei 



We now substitute expressions for A 12 i given by Equations (A31) and (A32) into 
Equations (A23) and (A24). Then, we replace the values Fji and $ iie i by their expressions 
through Ei which are given by 



1 jlx 



rrij 



rrii 



E± X + VjQy ( — 



E 



hi 



dE lx dE 



dy 



dx 



(A42) 



and 



*ii = - 



rrii 



9e 







d 



dx \dt 



d 



ViOy- 



d 



dy J \dt 



d \ 1 dE 



dy 



dx 



dE Xx dE lv \ q e d „ 
m P \ dy dx 



m e dt 



(A43) 
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When obtaining Equations (A42) and (A43), we have used Equations (A2) and (8). As a 
result, we will have the following equations for v ilx and vn y \ 
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For the electron velocity, using Equations (A26), (A27), (A34), and (A35), we obtain 
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and 
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A. 9. Perturbed plasma currents 

We now make use of obtained ion and electron velocities to find perturbed plasma 
currents j pilx = (qin i0 v ilx + q e n e0 v elx ) and j plly = (qin i0 v ily + qm^y + q e n e0 v ely ) in the 
general form. From Equations (A44) and (A47), we will have 
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and Upj = {Anrijoq 2 /rrij} is the plasma frequency. The values Aj je , /^, e , and are given 



by Equations (A46) and (A49). Using Equations (2), (A44), (A45), and (A48), we further 
find 
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We can rewrite Equations (A50) and (A52) in the form 
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where the components of the plasma dielectric permeability tensor are given by 
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Using Equations (A51) and (A53), we can find e p uj in specific cases. 
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B. Appendix 

B.l. Perturbed velocity of cosmic rays 

The linearized Equation (5) for the cold, nonrelativistic, T cr <C m cr c 2 , cosmic rays 
takes the form 

( d d \ f 2 Merger \ ^Pcrl -p-, Qcr /t~>i\ 

TcrO 777 + U crTT ^crl + 7cr0 7~ V crly = h F crl H V crl X B , (Bl) 
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When obtaining Equation (Bl), we have used that u cr is directed along the y-axis and 
7cri = lcro u crV C riy/c 2 , where 7 cr0 = (1 — u 2 r /c 2 ) . From Equation (Bl), we find the 
following equations for v CT \ x<lj : 
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where w ccr = q cr B Q /m cr c is the cyclotron frequency of the cosmic ray particles. Solutions of 
Equations (B3) and (B4) have the form 
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B.2. Equation for perturbed cosmic ray pressure 

From Equation (6) in the linear approximation, we obtain the perturbed cosmic ray 
pressure 
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Using the linearized continuity equation (2) for cosmic rays and expression for 7 crl , we find 
that p cr i is given by 
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From Equations (B5) and (B6), we obtain the expression for V • v cr i which is substituted 
together with the velocity v cr i y into Equation (B9). As a result, we have 
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where c scr = (p cr o^ cr / m cr n cro) 1 ^ 2 is the cosmic ray sound speed defined by the rest mass and 
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B.3. Equations for cosmic ray velocities via F cr i 



Substituting Equations (BIO) and (Bll) into Equations (B5) and (B6), we find 
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B.4. Equations for cosmic ray velocities via E x 



From Equation (B2), we have 
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Substituting Equation (B15) into Equations (B13) and (B14), we obtain 
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The operators Li crx ^ y , L 2cr , and L 3crxy containing in Equation (B18) are given by Equations 
(B7), (B 11), and (B12), respectively. 



B.5. Perturbed cosmic ray current 

We now find the components of the perturbed cosmic ray current j cr i x = q cr n crO v crix 
and j cr iy = q cr ncro v criy + q cr n cr iu cr . Using Equations (B16) and (B17) and continuity 
equation (2) in the linear approximation, we find 
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The components of the dielectric permeability tensor are the following: 
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